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Abstract. We give a multidimensional generalisation of the complete set of Bell-correlation inequalities given 
by Werner and Wolf in [26], and by Zukovirski and Brukner in [27], for the two-dimensional case. Our construc- 
tion applies for the n parties, two-observables case, where each observable is d-valued. The d'' inequalities 
obtained involve homogeneous polynomials. They define the facets of a polytope in a complex vector space of 
dimension d". We also show that these inequalities are violated by Quantum Mechanics. We exhibit examples 
in the three-dimensional case. 
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Tjj- ' 1. Introduction 

(n: 

The search for Bell inequalities has been the subject of a lot of work. Let us recall briefly what the matter is. 
^ Assume that a physical system is made of n subsystems. For each subsystem, a set of m different observables 

f^ I is considered. The outcomes of each of the nm observables belong to a set of cardinality d. The problem is 

to find inequalities which must be satisfied when a local realistic model is assumed. 
Ch ' The first such inequalities were provided by Bell [3] for the case {n,m,d) = (2,2,2). It was also shown 

2 ! that Quantum Mechanics violate these inequalities. The CHSH inequalities given in [5] were shown in [7] to 

be a complete set for the case (2,2,2). This means that these inequalities provide necessary and sufficient 

conditions for the existence of a local realistic model. 

The authors of [26] and of [27] gave a complete set of 2^ Bell inequalities for dichotomic observables, 
►^ . with arbitrary number of parties (case (n, 2,2)). The structure of these inequalities was further studied 

ly^ ' in [21], where a recursive method to compute Bell inequalities is also given. The tool for this construction 

lO ■ was the Walsh-Hadamard transform of Boolean functions. See also [24] which gives some insight and useful 

Cn ' details. 

(N 

l> 

O 



A method to obtain a complete set of dichotomic Bell inequalities was given in [19]. It has notably been 
t_^ I used to exhibit a complete set for the case (2, 3, 2). 

The multidimensional case has also been considered in numerous references. Reasons to explore beyond 
the two-dimensional case include that multidimensional entangled quantum states are known to be more 
resistant to noise, and that they can lead to stronger violations of local realism [13]. Also there are specific 
uses of the tridimensional case for quantum cryptography [14]. The pioneer work for multiple outcome 
^\ ' Bell inequalities was [6], where a family of multidimensional Bell inequalities, that generalize CHSH, was 

obtained. JVIoreover, these inequalities have been later proved tight [16]. 

However, no complete set has been given yet, beyond the two-dimensional case. 

Instead of the joint probabilities used by many authors for the multi- or three-dimensional case ([1], [6], 
[15], [19]), wc study the correlations between different observables using correlation functions. In general, if 
Xi{X) and Yj{X) are the values obtained by party i for the observable X and by party j for the observable Y, 
the corresponding correlation is given by 

' X^{X)Yj{X)p{X)dX 



A 



where A is the domain of the hidden variables A and p with J. p{X) dA = 1 is a density function. These 
correlation functions have been widely used for the study of the two-dimensional case, where the outcomes 
belong to {±1}. We use the same correlation functions also for the multidimensional case, but the outcomes 
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2. Multidimensional discrete Fourier transform 

arc now d-th roots of unity in C This approach has yet been considered for the d = 3 case in [9], [14], [15], 
[25]. 

We use a geometrical approach. Froissart [8] has apparently been the first to do so, and then the authors 
of [10] independently. It was shown in [18] that the local-realistic domain is a convex polytope (for joint 
probability distributions). The polytope corresponding to joint probabilities and the one corresponding to 
correlation functions are strongly related because of the relation E{X) — 2p{X = 1) — 1 between expectation 
values and probabilities, in the case d ~ 2. The polytope we consider belongs to a complex vector space of 
dimension d". 

Our inequalities are tight. This means that they define the facets of the polytope. The problem of 
obtaining all the (tight) inequalities was only solved in the two-dimensional setting ([7], [19] with joint 
probabilities, [26] and [27] with correlation functions). 

Our inequalities involve products and powers of observables, arranged in homogeneous polynomial ex- 
pressions. Powers of observables have already been used in [25]. It turns out that the method developed 
for (n, 2, 2) generalizes pretty well for the multidimensional, two-observables per party case, by means of 
multidimensional discrete Fourier transform. With this tool, we are able to give a complete set of tight Bell 
inequalities for the case (n, 2,d). 

In this paper, we first presents background about multidimensional Fourier transform (DFT for short). 
Then we recall some facts about the duality of polytopes in (finite dimensional) Hilbert spaces and study 
some useful relations between DFT and duality. Then we produce d"^ Bell inequalities which generalize 
those obtained in [26]. We study the polynomials involved in these inequalities and give some facts about 
the symmetries observed. Then we prove that our Bell inequalities form a complete set of tight ones. In 
section 7, we explain how violations of our Bell inequalities by Quantum Mechanics can be computed and 
observed. Finally we explore the case d = 3. 

2. Multidimensional discrete Fourier transform 

There arc numerous references for the discrete Fourier transform. One of them is [4]. However, we give here 
all the material we need for our purposes. 

Maps from Z^ to the set of d-th roots of 1 

The main tool for the classification of dichotomic Bell inequalities is the Walsh-Hadamard transform for 
Boolean functions. For our generalisation of the dichotomic case, we will use d-valued functions and multi- 
dimensional discrete Fourier transform. 

There are two equivalent ways to define Boolean functions: it can be a map F from {0, 1} to {0, 1} 
(additive convention), or a map / from {0, 1} to {1, —1} (multiplicative convention). The equivalence is of 
course given by / = (—1)-'^. The multiplicative convention is more comfortable when dealing with Walsh- 
Hadamard transforms. We also adopt a multiplicative convention, and the considered functions will take 
their values in the set 

11 = {l,w,...,cj'^"^} where w = exp(2z7r/rf). (1) 

We put Zd = {0, 1, . . . , d— 1} and denote by Z^ the set of n-tuples with components in Zd (rf, n E W). Also, 
we denote by J" or yd,n the set of maps from ZJJ to 11. There are d"^ such functions. 

The DFT 

Let / be a map from Z^ to the complex field C (or to IX as a particular case). The (multidimensional) 
discrete Fourier transform of / is the map DFT/ = /, also from ZJJ to C, defined by 

/(ri,...,r„)- ^ c.'-^^^+-+'^"*"/(si,...,s„) (2) 

si,...,s„GZtj 

or, written in compact form, /(r) — X^sPZ" ^^'^ f{s) where r ■ s — X^ILi ^«'** ^^ ^^"^ standard scalar product 

d 

of the n-tuplcs r and s. 



We denote as Hd the matrix {uj^^)o^i,j^d- 
D — d", given by 



2. Multidimensional discrete Fourier transform 

-I. The n-th tensor power of H^ is the D x D matrix, with 

The matrices i?®" can be built up from blocks using recursion on n: 

(1) and /fr = (^^-''^r"'W,,^„-i- (3) 



f7-®0 

^d 



These matrices are a generalization of the usual Hadamard matrices which are obtained in the special case 
d = 2 (hence lu = —1). 

A map / from Z^ to C can be identified to the vector of its values (/(s))sgz". The (column) vector of 
the values of / can be obtained applying the matrix i?®" to the (column) vector of the values of / : 



/ /(0,0,...,0) \ 
/(1,0,...,0) 



m 



( /(0,0,...,0) \ 
/(1,0,...,0) 

V/(d-l,...,(i-l)/ 



V/(d-i,...,d-i)y 

Hence, the map DFT: / h^ / is a linear map from C' to itself. 

Inverse DFT 

Let also define iJ^®" the matrix (w^''")r,isez"- It can be checked that 



TT*0n rr(g)n mj 



Hence, the inverse transform DFT is obtained by 



/(si,...,s„) = — J2 ' 

ri,...,r„eZd 



-(I'lSiH hr„s„) 



f{ri,---,Sn) 



or, in compact form, f{s) ^ ^ EreZ" ^ '""fir)- 

d 

In the particular case d = 2, the multidimensional discrete Fourier transform is the Walsh-Hadamard 
transform of Boolean functions: 

2;e{0,l}" 

(using the multiplicative convention: f(x) G {1,-1}). 

Some easy results 

Some easy results can be derived from the definition given by Equation (2), between the discrete Fourier 
transforms of two elements of 3^d,n which are related in some way: 

2.1. — Proposition. Put f = DFT/ and g = DFTg where f and g belong to 3^d,n- 

(a) Ifg{s) = fi^s) for all s e Z^, then g{r) = /(-r) for ah r e Z'J. 

(h) If g{s) = /(— s)* for all s e Z^, then g{r) — f{r)* for all r ^17^ (* denotes complex conjugation). 

(c) Let S e Z^. If g{s) ~ ,f{s + S) for all s G Z^ (addition in U^ is assumed component-wise and modulo d), 
then g{r) = uj-'"-^f{r) for ah r G Z^. 

(d) Let 5 e Z'J. Ifg{s) = oj^'^'fis) for ah s e ZJ^, then g{r) = f{r + S) for all r e Z^. 



3. Convex hulls 

(c) Let a be a permutation of the set {1, . . . , n}. For s = (si, . . . , s„) G Z^, we use the shorthand notation 
<j{s) = (s,(i), . . . , s,(„)). Ifg(s) = f{a{s)) for all s e Z^, then g{r) - f{a{r)) for ah r G Z^. 

Proof — We show only the last two assertions and leave the first three to the reader. Assume that 
g{s) = LU^-'fis) for all s G 17^. Then 



9 



(r) = Y. ^'^"•9(^) = E ^''■^^'"/(^) 



for all r e Z^\ Hence, 



9{r -5)^J2 ^^'^^'^■^^'■VG^) - E ^'^"/(^) ^ />)■ 



sei 



This proves assertion (d). Assume now that g{s) = f(^a{s)) for all s e Z^'. Then, for all r e Z^", 
= yj o;'"'^ '"'/(s) because cr induces a permutation on Z^. 



Hence 



3('^~'M) = E ^'^"^'"^■""''V(s) - E ^''■'/(^) - /w- 



seZJ S6Z2 

This proves assertion (e). D 

3. Convex hulls 

Let D e N. We denote (,5, 7) = J2i=iPiH the usual Hermitian inner product in C^. The complex 
vector space C^ can also be viewed as a vector space over M, with dimension 2D. Each element /3 € 
C'^ can be alternatively written as a D-uple (/3i, . . . ,/3_d) of coordinates belonging to C or as a 2D-uple 
(xi, 2/1, . . . , xdtVd) of coordinates belonging to M, with the relations j3k = Xk +iyk- Recall that the real part 
of the inner product (•, •) is nothing more than the usual scalar product in R^^: 

„ ,0 , o ^T" o* ^r , ^\ ii ^k = Xk + iyk i^nd -fk = Zk + itk, 

Re(/3,7) =Re > /3fe7fc = > (a;feZfc + yfeife) 

^ ^ with Xfc,2/fc,Zfc,tfe e M. 

Let S* be a subset of C''^. The convex hull of S is the set 

Hulls' := < y^ PkPk with (3k £ S and pj. e IR+ such that V^Pfe = 1 f- 

The dual (or polar) of the set S* is, by definition, the set 

T = 5° :={7eC^ |Re(/3,7) ^ l,V/3€5}. (4) 

When S* is a polytope containing 0, the vertices of the dual T correspond to the facets of S. To be precise, 
7 is a vertice of T if and only if the hyperplane defined by the equation Re(/3, 7) = 1 contains a facet of S. 
The following result holds (the bipolar Theorem, see [22]): 

3.1. — Theorem. For any subset S ofC^ containing 0, the dual S°° of the dual of S is the convex hull 
ofS. 



3. Convex hulls 



The hull of U and its dual 



We assume here d > 2. The convex hull of the set U is a regular polygon. The dual of U is also a regular 
polygon with d vertices (see Figure 1): 



3.2. 



Lemma. The dual U° of U (with d > 2) is the polygon with vertices set: 

1 



V = 



2k + 1 

■ exp I ; ITT 



fc = 0,.. 



1 



cos(7r/(i) \ d 

Proof — For I3k == exp (^) e U and 7, = exp f(21+l}iL^ 

Re ( exp{-2kiTT/d) exp {{21 + l)iTr/d)) cos ((2^ + 1 - 2k)TT/d) 



cos I ^ ) e V we have 



Re(/3fe,7/) 



cos(7r/(i) 



cos(7r/(i) 



Thus, Re(/3fc,7/) — 1 when k — I ot when k = I + 1 (the vertice ji of 11° corresponds to the edge Si — {/3i, 
/3;+i) of Hull II). For the other values of k, we have Ke{f3k,ji) < 1 because /3fc is in the half-plane delimited 
by di and containing 0. D 





Even example d ~ 4 Odd example d = 5 

Figure 1. The boundaries of the convex hull of U (solid) and its dual (dashed) 



3.3. — Lemma. DeGne p = exp(«7r/(i). For each (3 E Hull II, the following inequality holds: 

Re(p/3) :^ cos(7r/d). 

Hull II. Hence pf3 — cos{TT/d)j for some 7 e 11°. Thus, 



Proof — From Lemma 3.2, we have U° — — t-tjy 

' cos(7r/a) 

Re(p/3) = cos(7r/(i) Re(7). But we have Re(7) = Rc(l,7) ^ 1 because 1 e U. The result follows because 
cos(7r/(i) > (we assumed d > 2, note also that case d = 2 is trivially true). D 

Duality and DFT 

As in Section 2, we put D = d". The map DFT is linear and its matrix U = i/^" (in the canonical basis 
of C^) satisfies U'^U = DI, where U'' is the conjugate transpose of U. This has some useful consequences. 

3.4. — Lemma. Assume that /3,7 e C^, and put /3 = DFT/3 and 7 = DFT 7. We have (/3,7) = D{l3,"f). 

Proof — If we identify /3 and 7 with the column vectors of their coordinates in the canonical basis we can 
write: 

(/§, 7) = /3^ ^ {UpyUj ^ P^U^'U"/ ^ L>/3^7 = D{p, 7) 

as claimed. D 



4. Homogeneous Bell inequalities 

3.5. — Proposition. Let T be a polytope in C^ containing 0, and denote by T its image (which is also a 
polytope, by linearity of DFT) under the map DFT. We have the following relations between their duals: 

f° =Dr°. 

Proof — For /3 e C^, we have /3 e r° if and only if (/3,7) ^ 1 for all 7 G F. From Lemma 3.4, this is 
equivalent to {$,j) ^ D for all 7 e F. This condition can be written (-^/3,7) ^ 1, or therefore jj$ G F". 
Finally, it is equivalent to (3 £ DT°. D 

4. Homogeneous Bell inequalities 

Le n be the number of parties. For each party, we consider two observables, denoted by Ai and Bi (for 
1 ^ i ^ n). The outcomes of each measure are assumed to belong to the set U defined in (1), with d ^ 2. 
Recall, from the identity 

1-X'^ = (1- X){1 + X + X^ + ■■■+ x"^-^), 

that the roots of the polynomial 1 + X + ■ ■ ■ + X'^^^ are the elements of U \ {1}. Recall also that J2ui£U ^'^ 
evaluates to d when fc is a multiple of d but is zero otherwise. If 0^,6^ G U, there exists an integer r^ G Z^ 
such that Gi/bi ~ lu^\ Let also s^ G Z^. Then 

af-^ + Lo'^a^-% + ■■■+ uj^^-^^'^h^l-^ = a-^-^l + lo'^-"^^ +■■■+ Jd-i){s,-r,)) 

j af^^d if r, = s,, 

(^ otherwise. 

Let now / be any map from Z^ to U. We have 

n 

J2 /(^) n("^' + ^''<"'^» + • • • + ^'"'^^ a ■"'"'"' ^^r + • • • + LJ^''-^'>'^bt^) - ud" (5) 

where u G II, because in this sum, exactly one term is non-zero (the one corresponding to Si = Vi for each i). 
If we expand the products in (5), we get 

n 

t.d"=^/(.)5]n^""«'"'""^r 

n 

= J2 fir) a'' where / = DFT /. 

Now, if the ai and hi are random variables we can write, about expected values: 

Y^ f{r)E{a'^)ed'' HuhU. 

rezs 

From Lemma 3.3, we obtain: 

Re fp J2 /W^K)) «; rf"cos(7r/rf)- 

When d > 2, this also can be written: 

We call these relations homogeneous Bell inequalities. There are d^ of them. 



5. Homogeneous Bell polynomials 

5. Homogeneous Bell polynomials 

Wc now study the polynomials in 2n variables Ai and Bi (for I ^ i ^ n) which are involved in the 
homogeneous Bell inequalities. Some Bell polynomials where defined in [26] for d = 2. As a generalisation 
to the multidimensional case, we define the homogeneous Bell polynomials to be 

n 

Vf = J2 /»^'' ^^^'^ ^'' ■■- n ^'"'"'^•^r (7) 

where / is any map from Z^ to U. Let us denote ^d,n the set of these polynomials. Each element of O^d.n 
is a homogeneous polynomial of degree n{d — 1). Note that in view of Section 7, we consider IP/ as a non 
commutative polynomial. More precisely, each Ai is not assumed to commute with Bi, while Ai and Bi do 
commute with Aj and Bj for i ^ j. 

As in [22] where the case d = 2 is handled, we give a recursive construction of the homogeneous Bell 
polynomials. This construction is a direct consequence of Equation (3). If 5*0, ... , Td-i sue homogeneous 
Bell polynomials in the 2(n — 1) the variables Ai, Bi with 1 ^ i ^ n — 1, then we get a homogeneous Bell 
polynomial in 2n variables by the d-ary operation c><i: 

r„=0 ^ t=0 ^ 

Conversely, every element of the set ^d,n can be obtained this way. 

For example, with d = 2, the polynomials obtained are ±1 for n = 0, ±2Ai and ±2Bi for n = 1, and 
±AAiA2, ±2{-AiA2 

±4AiB2, ±2(^1^2 

±451^2, ±2(^1^2 

±4BiB2, ±2{AiA2 

for n — 2 (we recognize the polynomials involved in the CHSH inequalities). Examples for d = 3 will be 
given in Section 8. 

Symmetries 

The set Jid.n of homogeneous Bell polynomials has some symmetries we briefiy discuss now. They are 
consequences of Proposition 2.1. 

a. If the maps / and 5 G Z^ are the same, up to the order of their arguments: 

g{si,...,Sn) =/(s<t(i),...,Sct(„)) for all s e Z^, 
for some permutation cr, then the polynomial Tg can be obtained from Vf by changing each variable Ai 
(resp. Bi) to A^r^) (rcsp. B^ii\). This symmetry corresponds to the fact that the n subsystems arc indistin- 
guishable. 

b. If, for some iq, 

g{s) ^ uj-'-o f(s) for all s G Z^, 

then, from Proposition 2.1, we have g{r) — f{r — S) for all r E Z^, where 6 — (0, . . . , 0, 1, 0, . . . , 0) has its 
only non- null component at index iq. Hence, we obtain 

T,= J2 9ir)A^ = E /(^ - '^)^'" = E fi^W^'- 

This shows that we obtain Tg from Vf by the circular monomial substitution 

Af-^ — ^ Af-^B,„ -^ > A,„Bf-^ -^ Bf-^ -^ Af-\ 

Also, the set IH^.n is invariant, under the swap operation Ai^ o Bi^ (this can be algebraically checked with 
the help of Proposition 2.1(a)). Hence, for each iq, the set 'Kd.n is invariant under the action of the dihedral 
group of order 2d over the monomials made of the variables A^^ and Bi^ . 

c. Of course, the set ^d,n is also invariant under multiplication by cj, and by complex conjugation 
(Proposition 2.1(b) can be used to check this latter fact). 
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6. The classical domain 

6. The classical domain 

Wc now show that the homogeneous Bell incquahties obtained in Section 4 are tight and completely character- 
ize a local realistic model, for n € N* parties, m — 2 observables for each site, and d-outcomes measurements 
with d> 2. 

The values a^ and 6^, when a local realistic model is applied, of these 2n observables are assumed to 
belong to the set U. We consider the monomials 

n 
i=l 

which appear in homogeneous Bell polynomials. There are D ~ dJ^ of them. For each experiment, the data 
set of the values obtained for these monomials form a vector ^ = (a'')sgz" in C^. Our aim is to show that 
the domain accessible to the expected values of ^ is the polytope defined by the inequalities (6). 

The polytope 17 

Put 

i^r = K-'')seZ5 e C^ for each r e Z^. 

The (i"+^ = dD vectors u^n for u G U and r G Z^' are all distinct. In a local realistic model, each experimental 
data set assigns a value 



n«' ""^r-^«^'^-""=^ 



at^uj"--' 



to each monomial vl" where w'"' — hi/ai (for 1 ^ i ^ n). Thus, the vector ^ obtained from experimental data 
is one of the vectors u^n where u — Y\^^i cl^~ , and r — (ri)i^i^„ with the r^ just defined. 

Conversely, it is possible to design classical experiments which assign independently any value in XL to the 
2n variables and which assign any u^^ to the data set vector ^. Then, if the values assigned to the variables 
follow some probability distributions, expected values for the vectors ^ obtained, are convex combinations 
of the u^r- Hence the classically accessible region for ^ is the convex hull of the u^^, which will be denoted 
by Q, as it was in [26] for the case d = 2. The domain il is a polytope in C^ and has dD vertices. Notice 
that il has a d-order symmetry: loQ. = i7. 

The polytope H = DFT"^ Q 

We can find all the inequalities defining the facets of the polytope J7. They will be the d^ homogeneous Bell 
inequalities (6) we obtained in Section 4. 

Let {'Ks)sei"' be the canonical basis of the complex vector space C^. The discrete Fourier transform 
maps the tt^ to the ^g. We consider the following polytope: 

n := HullJMTrs I u e 11, s e Z^}. 

Then Q, is II, the image of 11 under DFT. To find the facets of 51, we have to study its dual. But from 
Proposition 3.5, 

n°=II° = ^IP. (8) 

Let's first study 11°. 
6.1. — Proposition. The vertices of the polytope 11° are the /3 = {/3i, . . . , /Sg) such that 

Ps — 7 — j-irfis) where f is any element ofJdn- 

cos[n/d) 

Proof — By definition, 

n° ={/3eC^ |Re(/3,u7r,) s$ 1, Vm G U, s e Z^}. 



7. Violations by Quantum Mechanics 
Using the rf-ordcr symmetry ofll", and using [jS^wKs] — u{/3,t:s), we can write 

We are interested with the extremal points of 11° . These are obtained when {(3, tTs) are in a corner of 11° (see 
Lemma 3.2): 



COS(TT/d) 



Hence, there exists / € 3^d,n such that: 



/3: = (/3,7r,) = ^--/(s) forahseZS. 
cos(7r/aj 

But n° is symmetric under complex conjugation. Hence we can change l3* for /3s. 
The dual of 17 

6.2. — Theorem. The vertices of the polytope il" are given by 

-(/>)) ^^„„ for/eJ,,„. 



d"cos(7r/d) ^^ ^^rez^ 

Proof — The result follows from Equation (8) and Proposition 6.1. D 

To end this section, note that the inequalities (6) can be written 

Re(/3„0O with /3;- = ^;r^4^(/»).,.2 -^ ^-(^K)).ez,- 

Hence the theorem just obtained shows that our homogeneous Bell inequalities define the facets of the 
polytope V,. Thus they form a complete set of tight Bell inequalities. 

7. Violations by Quantum Mechanics 

At this point, we have only considered local-realistic models. The polytope 51 we have made explicit using 
homogeneous Bell inequalities is the domain accessible with such models. However, the primary aim of Bell 
inequalities was (at least historically) to compare local-realistic theories with Quantum Mechanics. The main 
success of the original and CHSH Bell inequalities, was due to the fact that Quantum Mechanics violate 
them, hence they provided the proof that quantum indeterminacy cannot be explained by hidden variables. 
We now show that Quantum Mechanics also violates homogeneous Bell inequalities, and that this fact could 
be, in principle, checked by experiment. 

There exists a difficulty in our setting, which did not appear in the d — 2 case. Multidimensional 
homogeneous Bell polynomials involve products of variables, some of them corresponding to observables of 
the same party. In Quantum Mechanics, such observables corresponds to non commuting operators and 
their values cannot be simultaneously obtained, and this prevents to observe violations this way. However, 
there are important cases where such products of observables are themselves observables. This is our key 
tool now. 

Generalized Pauli matrices 

We use the following multidimensional generalization (found for example in [23] and [11]) of Pauli (or spin) 

matrices. Let 

/O ••• 1\ /I ••• \ 

' w_i 

cj •• 



X 



1 ■■ ■■ : 

: '•. '•. 
VO ••• 1 Oj 



d-i 



yj ji + 1 mod d)(i| and 



d-l 



Y,o^'\^){^\ 



\0 ••• w'^-V 

where the kets \0), . . . ,\d — 1) form an orthonormal basis of C* (in fact, the eigenbasis of Z). The matrices X 
and Z have order d and satisfy ZX = ujXZ. The generalized Pauli matrices are the following d+1 unitary 
matrices: 

Z, X, XZ, ..., XZ'^-\ (9) 

The following two results are easy to show. The first one is about eigenvalues as these are the possible 
outcomes of measurements in Quantum Mechanics. 



7. Violations by Quantum Mechanics 

7.1. — Proposition. Let k be an integer. The eigenvalues of XZ^ are the lo^ (with ^ j ^ rf — IJ when 
d is odd or when k is even. They are the puj^ , with p — cxp^in / d) , when d is even and k is odd. 

Proof — By expanding the characteristic polynomial of XZ'^ along the last column, we obtain : 
detiXZ'' - XI) = (_i)d-i^fe(d-i)^fc(o+i+-+(d-2)) _ xi-\^-^) 

= {-\f - (^^l-^'i^k(d-l)+k{d-l){d-2)/2 ^ ^_^-jd _ ^_^-jd^kid-l)d/2 

( ^xd _^ ^k(d-i)/2 ^i_xd when d is odd, 

I ^d _ ^_-^)fc(d-i) ^x"^ - {-if when d is even. 

Hence the eigenvalues are the solutions of equation A'' = 1 when d is odd or k is even ; and of A'' = — 1 when 
d is even and k is odd. D 

7.2. — Lemma. For any integers k, e, the following relation holds 

Proof — We leave it to the reader. It can be done by induction over e, using ZX = ljXZ. D 

Unitary observables 

It is shown in [2], and also in [20], that when d is a power of a prime, the bases consisting of the normalized 
eigenvectors of the rf + 1 Pauli matrices given by (9) form d+1 Mutually Unbiased Bases [12]. Paterek [17] 
explains that these generalized Pauli matrices can be used as unitary observables, instead of more classical 
Hermitian operators. (Note that for the usual case d — 2, the matrices X and Z are Hermitian as well as 
unitary.) These are clearly the operators we need, as we considered complex valued observables. 

To determine quantum violations of homogeneous Bell inequalities, we have to evaluate expected values 
of operators of the form X'^^^^^Z^, for ^ r ^ d — 1. These are unitary observables, and it should be 
possible to directly obtain outcomes of them (without measuring outcomes of X and Z). At least, we can 
rely on the better known generalized Pauli operators using the following Proposition. 

7.3. — Proposition. Let ^ r ^ d—1 and assume d prime. It is possible to experimentally obtain values 
for x'^^^^'' Z"^' in order to compute the corresponding expected value. 

Proof — When r = d — 1, just make a measurement with operator Z on each sample, and raise the 
outcomes to power d~ 1. We now assume that ^ r ^ d — 2. Thus l^d— 1 — r^d— 1. Asdis prime, 
then (d — 1 — r) is invertible modulo d and it is possible to find an integer k such that 1 ^ fc ^ d — 1 and 
k{d — 1 — r) = r modulo d. From Lemma 7.2, we get 

I -^7- r,k\d-l-r _ k(d-l-r){d-2-r)/2 vd-l-r yk(d-l-r) _ k(r+l){r+2) /2 v-d-l-r yr 

Hence, we have to make a measurement with operator X Z'^ on each sample, raise the outcomes to the power 
d — 1 — r, and multiply the results with a;^'=('"+i)(''+2)/2, □ 

Now, we are able to compute some violations of homogeneous Bell inequalities by Quantum Mechanics, 
with the quantum operators X and Z in place of the classical operators Ai and Bi respectively. Hence, we 
consider the following quantum counterparts of our homogeneous Bell polynomials (7): 

n 

Qf = Y. ■hr)AQM where A^q^ := (S)iX'-'-'''Z^'). 
A quantum state \(j)) will violate the corresponding homogeneous Bell inequality if the condition 

is not satisfied. Our short study of the case d = 3 will indeed exhibit cases where such violations occur. 
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8. The case d = 3 

8. The case d = 3 

We illustrate our results with the first multidimensional case: d — 3 (sometimes called trichotomic). Note 
that the factor l/cos(7r/(i) in Equation (6) is maximal in this case, and this might lead to higher violations. 

DFT 

Here, w = exp(2i7r/3) and 
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The huh of XL and its dual 

The hull of U is the triangle with vertices 1, w, uP' and its edges arc defined by the three inequalities 

x + \/3y<l, -2a; < 1, x-%/32/:^l. 

Hence, the dual 11° has vertices 1 + i\/3, —2, 1 — iVS, which are obtained from the vertices of HuUU by 
multiplication by exp(i7r/3)/cos(7r/3) = — 2cj2. 

Beh polynomials 

We did some computations, with the help of the Magma computer algebra system. For the (virtual) case 
n = 0, the trichotomic Bell polynomials are the constant polynomials 1, lo and lo^ . For n — 1, there are yet 
27 homogeneous trichotomic Bell polynomials. Instead of listing them all, we give for them the following 
compact expression: 

u{3M + (v - 1){A'^ + AB + B^)) (10) 

where u,v e {1,lo,lo'^} and M e {A'^,AB,B'^}. 

For n — 2, there are 19683 homogeneous trichotomic Bell polynomials. Among them, 18792 are ir- 
reducible polynomials. The number of elements in ^{3,2 with only real coefficients is 81 (the aim of this 
criterion here is just to reduce the list size). We can list them, up to the symmetries discussed in Section 5, 
as there remain only 4 ones: 



3{AlAl-AlBi 



9AlAl 
2A1B1A2B2 + AiBiBJ - BJAl 



3{-AlA2B2+AlBl 
3{2AlAl - AIA2B2 - AJBi 



AiBiAi + 2AiBiBi 



BfAJ 



A1B1AI+A1B1A2B2 



BIA2B2) 
^BlA2B2) 
^AiBiBl). 



We found also that there are 243 elements in ^{3,2 up to these symmetries. 

Bell inequalities 

The factor ^„ coMt^ Id) ^i^h appear in Inequalities (6) is in this case — 2a;2/3". By changing / to w/, we can 
remove the w^ to obtain the following homogeneous trichotomic Bell inequalities: 



Re ( T^ E /('^)^(«'') ) < 1 for each / e J; 



3,n- 



(11) 
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9. Conclusion 
Violations 

Yet the case n = 1 is especially interesting. The 27 homogeneous Bell polynomials fall in 3 classes according 
to the value of v in formula (10). The most interesting class is the one obtained with v — lj'^. In that case, 
the eigenvalues of the operator obtained are — 3^, — SCw and — SCw^ where ( = exp(2i7r/9). They do not 
belong to Hull II. In particular — 2Re(— 3C)/3 ~ 1.53209 > 1 and one can expect violations. This is indeed 
the case: consider the map / such that /(O) = lu and /(I) = /(2) = cj^. Then we have /(O) — ui'^ — 1, 
/(I) = /(2) = cj — cj^ and Equation (11) reads 



-- Re ({lu^ - l)E{a^) + (cj - iu^){E{ab) + E{b^))\ «: 1. 
The corresponding operator is 

oj-oj^ l-oj a;2 - 1 
cj2-l a;2-l u;2-l^ 

and has A — ~3( as an eigenvalue. Then we can find states \tp) such that — 2Re(i/'|(5j^|i/')/3 exceeds 1. For 
example, the state (|0) + 2|1) + 3|2))/a/14 achieves a violation of 19/14 ~ 1.357 and the state 



((44 + 50w)|0) + (76 + 9u;)|1) + (143 + 17cj)|2))/V25716 

achieves a violation of 1.53208. Non-locality is not needed to violate homogeneous Bell inequalities! Of 
course, this situation did not appear in dimension rf == 2, as there were only trivial Bell polynomials with 
n = 1. 

For n ~ 2, the best violation is obtained in 27 cases for which eigenvalues are 9(1 — w), 9(a;^ — 1), 
9{uj — u;^) and (with multiplicity 6). One of these cases is the following 

Q/ =3({uj^ -1)X^(^XZ+ (cj2 - l)XZ (^ X^ + {1 - u})Z^ (^ Z^] 

obtained from (11) with the map / whose vector of values is (u;^, (x',w^,a;,a;, l,^^, 1, 1). We expect violations 
of -2Re (9(a;2 - l))/9 = 3. Such violation is obtained with the state (|01) + |10) + uj\22))/V3. 

9. Conclusion 

In this paper, we defined homogeneous Bell inequalities and we showed that they correspond to the boundaries 
of the domain accessible with local-realistic models, for the general multipartite and multidimensional case 
with two observables per party. We studied homogeneous Bell polynomials and their symmetries. It turns 
out that the classical domain is the image under DFT of a polytope obtained from the canonical basis, and 
we used this fact to compute its dual. With this, we were able to show that the homogeneous Bell inequalities 
form a complete set. 

Then we considered violations by Quantum Mechanics, using the observables provided by generalized 
Pauli matrices. We showed that violations indeed occur, and exhibit some of them in the trichotomic case. 

The complex valued correlation function we used is a natural mathematical generalisation of the two- 
dimensional one. Fu in [9] argued that it has also a physical meaning, at least in Quantum Mechanics. It was 
a crucial and fruitful ingredient in the present work, and this raises interrogations about the precise extent 
of this physical meaning. Also, complex valued observables provided by the generalized Pauli matrices were 
a key tool for computing violations. 
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